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AIRBORNE LASER-RADAR STUDIES OF THE LOWER ATMOSPHERE 


By Samuel H. MeLfi and Joseph W. Stickle 
Langley Research Center 

SUMMARY 

In the present investigation a theoretical model of ruby-laser radiation backscatter 
with a wavelength of 0.6943 pm in the lower atmosphere has been developed. This model 
accounts for scattering from air molecules and atmospheric aerosols. Molecular scat- 
tering appears to predominate over aerosol scattering in the altitude range from 4 to 
14 km. Above 14 km, an aerosol layer increases the scattering by approximately a factor 
of 2 above Rayleigh (molecular) predictions. An airborne laser system has been devel- 
oped and flown from near sea level to an altitude of 11 km. The results of these airborne 
measurements agreed well with the model predictions. 

This investigation was initiated to investigate the possible correlation between 
molecular or aerosol density fluctuations due to clear-air turbulence and laser-radar 
return from the atmosphere; however, no correlation was obtained. This result, however, 
does not preclude using laser backscatter systems as atmospheric probes, since the sys- 
tem used in the present investigation was extremely insensitive to the detection of molec- 
ular or aerosol density fluctuations. By using the results of this investigation and the 
theoretical model, an analysis of a general airborne laser system has been performed. 
This analysis has indicated that use of an improved system with better sensitivity and 
which is within the state of the art would make the airborne laser system a valuable tool 
for investigation of atmospheric phenomena associated with molecular and aerosol density 
changes . 


INTRODUCTION 

Many investigations have been made on the scattering of light in the atmosphere. 

The first proposal for obtaining density profiles by measuring scattered light from a 
ground-based searchlight was made in 1930. (See ref. 1.) In 1937 (ref. 2), 1939 (ref. 3), 
and 1954 (ref. 4), molecular densities above 10 km were measured by using the technique 
of reference 1. More recently, with the advent of the high-powered ruby laser, light- 
scatter probing of the atmosphere has been extended to higher altitudes. It was reported 
in 1963 that atmospheric backscatter from a ruby laser had been measured up to an alti- 
tude of 140 km. (See ref. 5.) Subsequent investigations in which laser probing techniques 


were used are reported in references 6 to 10. In reference 9 "ghost" returns from 
apparently clear air when using a ground-based ruby laser were observed. In refer- 
ence 10 these ghost returns were attributed to scattering by layers of atmospheric aero- 
sols. By using an airborne laser, an inconclusive correlation between laser returns and 
clear-air turbulence was found. (See ref. 11.) 

Based primarily on the discoveries reported in references 9 and 10, a program was 
initiated at Langley Research Center to probe the atmosphere by using an airborne ruby 
laser. To evaluate the capability of a laser system to detect molecular and aerosol den- 
sity gradients, a theoretical model of light scattering in the atmosphere is developed 
herein. The model is then verified by comparing it with experimental data obtained while 
flying a laser -equipped airplane at various altitudes from sea level to 11 km. With this 
verification, the model is used to determine the capabilities of laser systems in general 
and to predict the performance of an improved system that is well within the present state 
of the art. 
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SYMBOLS 

effective receiver area, meters^ 

transmitted beam cross-section area at location of scattering volume, meters^ 
change in aerosol density, percent 
receiver bandwidth, hertz 

constants of aerosol distribution at altitude h 

velocity of light, 3.00 x 10® meters-second - ! 
depolarization factor 

total energy transmitted by laser pulse, joules 
anisotropy factor 

scattering function, steradian~l-meter“l 

Mie scattering function with 8 = 180° and A = 0.6943 micrometer, 
steradian - 1-meter _ 1 
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fM(S,^jh) Mie scattering function, s ter adian"l- meter “1 

fp^(h) Rayleigh scattering function with 9 = 180° and X = 0.6943 micrometer, 

steradian" 1 -meter" 1 

f R (0,X,h) Rayleigh scattering function, steradian - 1-meter" 1 

ff(h) total scattering function with 6 = 180° and X = 0.6943 micrometer, 

steradian" 1-meter" 1 

G gain of the photomultiplier 

g nondimensional size parameter of atmospheric aerosols, 2irr/\ 

h altitude of airplane, meters 

ij(e,X,L,p)l Mie scattering amplitudes of scattered light perpendicular and parallel, 

i 2 (0,X,t,p ) J respectively, to plane through direction of propagation of incident and 

scattered beams 

Ja,D anode dark current, amperes 

K parametric variable 

k wave number, 2ti/\, micrometer"! 

m change in molecular density, percent 

N e g number of photoelectrons produced at photocathode of receiver by sky 

brightness in a counting time t c , photoelectrons 

N e d number of thermal electrons (dark current) produced at photocathode of 

receiver in a counting time t c , electrons 

N e R number of photoelectrons produced at photocathode of receiver by Rayleigh 

scattering in a counting time t c , photoelectrons 

N e> t number of photoelectrons produced at photocathode of receiver by Rayleigh 

and Mie scattering in a counting time tc, photoelectrons 
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number of Mie scattering centers in scattering volume 
number of Rayleigh scattering centers in scattering volume 
number density, meter -3 

number density of Mie scattering centers at scattering volume, meter~3 

number density of Mie scattering centers with aerosol radius r at alti- 
tude h, meter“3 

number density of Rayleigh scattering centers at scattering volume, meter~3 
power incident on receiver, watts 
average power transmitted by laser, watts 
shape of aerosol particles 

receiver subsystem quality factor, rjyA Y jAB, second-meter^ 
range of the scattering volume in front of airplane, meters 
radius of atmospheric aerosols, micrometers 
signal-to-noise ratio 

signal-to-noise ratio due to a change in aerosol density 

signal-to-noise ratio due to a change in molecular density 

sensitivity of receiver subsystem, volt-microwatt - * 

counting time of receiver, seconds 

elapsed time as measured by receiver, microseconds 

voltage at oscilloscope, volts 

scattering volume, meters^ 
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pulse width, meters 

average molecular polarizability of air up to h = 80 kilometers, meters^ 
extinction coefficient due to Mie scattering centers, meter - 1 
extinction coefficient due to ozone molecules, meter “1 
extinction coefficient due to Rayleigh scattering centers, meter-1 
total extinction coefficient, meter - ! 
receiver optical efficiency 
field of view of receiver, mrad 
beamwidth of laser transmitter, mrad 

conversion factor from energy to photons at X = 0.6943 micrometer, 

3.51 X 10l8 photons-joule - ! 

receiver photocathode quantum efficiency, photoelectrons/photon 

scattering angle (measured between forward direction of incident beam and 
direction of observation), degrees or radians 

average index of refraction of the atmospheric aerosol 

curve constant of the atmospheric return, volts-second^ 

wavelength of incident light, micrometers 

optical bandwidth of the receiver, micrometers 

standard deviation 

differential cross section, meters^-steradian - ! 

Mie extinction cross section, meters^ 
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Subscript: 


differential Mie scattering cross section, meters^-steradian"! 
differential Rayleigh scattering cross section, meters^-steradian - ! 
width measured at half-power points of laser pulse, seconds 
polar coordinate, degrees or radians 
solid angle subtended by the receiver, steradians 


C 


correction for attenuation effects in the atmosphere 


THEORY 


A diagram of the general airborne backs cattering problem is shown in figure 1. 
Scattering of light in the atmosphere is produced by scattering by air molecules and 
atmospheric aerosols. Equation (1) gives a relationship between the power transmitted 
by the laser and the power incident on the receiver and is derived in appendix A: 

P r = — ^^ R (e,X)n R (h) + a M (e,A)n M (h)]e' 2 ^ t(h)R (1) 

2R 2 

There exists a general relationship between the scattering function, the differential cross 
section, and the number density which is given by 

f(e,A) = a(0,A)n (2) 


By using this relationship, equation (1) reduces to 


Pr- 


cE T A r[- f R (0jMi ) + fM ( 0jAjh) ] e - 2 ft (h)R 


2R 


(3) 


Equation (3) indicates that scattering in the atmosphere can be analyzed by considering 
separately Rayleigh scattering (including molecules and aerosols with diameters much 
smaller than the incident wavelength) and Mie scattering (aerosols with diameters on the 
order of the incident wavelengths). 
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The scattering functions for Rayleigh scattering fp,(h) and Mie scattering f M (h) 
have been derived in appendixes B and C, respectively, with 6 = 180° and 
A = 0.6943 jum and by assuming a polydisperse collection of spherical aerosols as 


fft(h) = (2.10 x 10“ 32 )n R (h) sr -1 m -1 

r 

f M (h) = ^6.0 x 10“ 9 jn^i(h) sr -1 m-l^ 


(4) 


These quantities are shown plotted as a function of altitude in figure 2. The molecular 
number density of the Rayleigh scattering function is a function of altitude, which can be 
obtained from reference 12. The aerosol number density of the Mie scattering function 
is developed in appendix D and given in table I. Atmospheric extinction coefficient /3^_ (h) 
is presented in appendix E. The solid curve in figure 2 is the sum of the Rayleigh and Mie 
scattering functions. It represents the total scattering function for clear air as a function 
of altitude. This model predicts that backscatter light at a wavelength of 0.6943 pm in the 
altitude range from 4 to 14 km is primarily due to air molecules. Above 14 km, an aerosol 
layer increases the scattering by approximately a factor of 2 above Rayleigh predictions. 


EXPERIMENTAL VERIFICATION OF THE MODEL 


Description of the Airborne Laser System 

A block diagram of the airborne laser system used to verify the theoretical model 
is shown in figure 3. This system can be conveniently divided into three subsystems: 
the transmitter subsystem, the receiver subsystem, and the display and recording subsys- 
tem. The location of the three subsystems is indicated on a photograph of the jet airplane 
used in this investigation. (See fig. 4.) Figures 5 and 6 are photographs of the equipment 
removed from the airplane. 

The transmitter subsystem consists of a Q-switched ruby laser, collimating optics, 
and an output monitor. The characteristics of this subsystem are summarized in table H. 
Temperature control of the laser head by a closed-loop cooling system maintains the 
laser wavelength at 0.6943 pm. The output monitor detects the laser radiation as it 
passes through the collimating optics. The output of the monitor photodiode is used to 
trigger the oscilloscope in the display and recording subsystem and is also electrically 
integrated to provide a measure of the energy in each laser pulse. 

The receiver subsystem consists of receiving optics, a temperature-controlled 
interference filter, and a photomultiplier tube. The characteristics of the receiver sub- 
system are summarized in table II. The display and recording subsystem consists of an 
oscilloscope and 35-millimeter camera. The oscilloscope has a bandwidth of 50 MHz, 
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which is the limiting bandwidth of the system, with a two-channel input. The camera has 
an automatic film advance which is capable of framing at 1 frame/sec. 

Calibration and Alinement 

The transmitter and receiver were calibrated before and after each flight, and their 
optical axes were alined periodically during the series of experiments. 

The laser transmitter was calibrated by using a thermopile calorimeter. Prior to 
the 28th flight test, the output power used in reducing the data was assumed to be an aver- 
age of the pre- and post-flight calorimeter measurements. Beginning with flight 28, a 
photodiode was utilized to record the output of individual pulses during a flight. The pho- 
todiode was calibrated by firing the laser into the calorimeter a number of times (10 
to 20) and recording the output of the calorimeter along with the integrated output of the 
photodiode. These data were plotted. A straight-line fit to the data was made which gave 
a constant relating total energy output of the laser to maximum deflection on the oscillo- 
scope of the integrated output of the photodiode. This output of the photodiode was 
recorded with each backs catter measurement; thus, the record gives both the scattered 
return and the output energy for each shot. 

The receiver sensitivity was measured by using a standard lamp, which was cali- 
brated by the National Bureau of Standards and was found to be 0.94 volts/pW. In order 
to perform this measurement, it was necessary to know the solid angle subtended by the 
receiver and its spectral bandwidth. With a measured half-angle field of view of 3 mrad, 
it was determined that the solid angle was 28 psr. The spectral bandwidth of the system 
was limited by the narrow-band interference filter. An analysis of the interference filter 
determined that the energy received from the standard lamp was predominately within the 
0.001175 fim band pass of the filter. 

The alinement of the system was performed by adjusting the direction of the laser 
beam until a maximum return was obtained from atmospheric scattering. 

Analysis of Data 

The data were analyzed and compared with results from the theoretical analysis. 

The airplane was flown at predetermined altitudes and the laser was fired at 1 pulse/sec. 
Figure 7 is a typical trace obtained. This probing was taken with the airplane flying at an 
altitude of 0.31 km over Williamsburg, Virginia (lat. 37°17' N, long. 76°45' W). Labeled 
on the trace is the peak of the integrated output of the photodiode monitor, crossover of 
the transmitter beam and the field of view of the detector, and the normal R - ^ falloff of 
the atmospheric return. 
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At each altitude, at least five probings were made. The output energy and the sig- 
nal return at -i-jiisec intervals for each probing were later read and recorded. These 

LX 

data were then analyzed in the following manner: Equation (3) can be written as 


Pr= ^r ft(h )e- 2 WR 


(5) 


2R^ 


Since for this system the range R is always very small, ~ 1; thus, 

f t (h) 


cEfjiA p 


r 2R 2 


( 6 ) 


Now, the signal voltage at the oscilloscope can be given as 

V(t E ) = P r x s r x 10 6 volts 


where 


R = M ^ x io 2 m 


By combining equations (6), (7), and (8) 


^ cE T A r s r 




f t (h) x 10 2 volts 


The general form of the V^t^j distributions can be written 
Combining equations (9) and (10) and solving for f^(h) yields 


(7) 


( 8 ) 


(9) 


( 10 ) 


f t (h) = | —4 x 10" 2 sr _1 m“l 

2 cE T A r s r 


(lla) 


With A r = 9.48 x 10“ 3 m, 


f t (h) = 


1.58k 

E'j'Sj. 


x 10" 8 


sr“lm"l 


(11b) 
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Thus, to arrive at the total experimental scattering function, it is necessary to measure 
the curve constant k, the total energy Eq>, and the receiver sensitivity s r . 

Error Analysis 

Two types of error are considered to influence the results of the experimental data. 
The first is an absolute error associated with the calibration equipment and technique, 
and the second is a random error arising from system fluctuations and data reduction. 
The absolute error is estimated to be ±25 percent based on ±10 percent error in the 
transmitter calibration due to calorimeter uncertainty, and ±15 percent in the receiver 
calibration. The random errors considered are estimated as follows: 


Root-mean-square 
error, percent 


Photomultiplier electron emission (shot noise) ±5 

Output monitor ±3 

Frequency stability of 400 Hz power supply ±2 

Reading error ±3 


These errors are considered to be independent of each other and, therefore, combine as 
the square root of the sum of the squares to give a 1-a value of ±6.9 percent. With a min- 
imum of 5 probings used to compute the average scattering function for each altitude, the 
overall random root-mean-square error reduces to <?/{$ (3.06 percent). The overall 
3-cr error then is 9.18 percent. Therefore, the values of scattering function presented in 
this report are considered to be absolute within ±34 percent. This percentage applies 
only to data obtained in flights 28, 29, and 31 for which the output power was monitored in 
flight. 


Results 

Figure 8 presents a comparison between the experimental results obtained on 
March 9, 1967, over Williamsburg, Virginia (lat. 37°17 T N, long. 76°45' W), and the cal- 
culated theoretical model. Each datum point represents the average scattering function 
for that altitude taken from five consecutive laser probings. For each probing the total 
transmitted energy was measured, and the curve constant k was obtained by computing 
the best fit t" ^ curve to the signal-versus-time data. (See fig. 7.) The experimental 
data are seen to be in very good agreement with the theoretical predictions. 

Figure 9 presents a composite of additional data taken from March 1966 to April 
1967 near Williamsburg, Virginia and over the local coastal region of the Atlantic Ocean. 
These data tend to support the theoretical model; however, there are two important points 
to be made. The first is that a measure of the output energy for each pulse is almost 
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essential for probings of this nature. For the data shown for all flights before flight 
number 28, the output energy was taken to be the average of measurements made on the 
ground before and after each flight. From flight 28 on, the energy was measured for each 
pulse while in flight. The results show that the scattering functions measured in flights 
28 through 31 are in closer agreement with the theoretical model than results from 
earlier flights. 

The second point is in regard to defining scattering of a clear atmosphere. Since 
the only basis for classifying clear air was from pilot's and observer's notes recorded in 
the flight log, there is obviously room for questioning the term "clear air." The three 
data points shown for flight 31 at an altitude of about 9 km were taken in a very light haze 
layer. While there is considerable variation of scattering function between the samples, 
all three are considerably greater than the data taken in the air judged to be "clear." 

Although the system was sufficiently sensitive to provide a measure of the atmo- 
spheric scattering function, it did not have sufficient performance capability to probe 
small horizontal scattering changes that might be associated with such phenomena as wind 
shear layers, temperature gradients, and possibly clear-air turbulence. 

APPLICATION OF THE MODEL TO LASER SYSTEM DESIGN 

Introduction 

In the preceding sections a model was developed to describe the backscatter of 
ruby-laser light in a clear atmosphere. In the following section the model has been used 
to generate a series of parametric curves to aid investigators in the evaluation of detec- 
tion capabilities for existing systems and to aid in the design of a system based on some 
predetermined detection criteria. The signal detectability of a system is discussed for 
three hypothetic atmospheric conditions: 

Condition I is total signal detectability; the molecular and aerosol densities are 
horizontally uniform and adhere to the model shown for the total scattering function. 

Condition H is molecular signal detectability; the molecular and aerosol densities 
conform to the model out to a horizontal range R; at and beyond R, an increase or 
decrease in the molecular density exists. 

Condition IH is aerosol signal detectability; the molecular and aerosol densities 
conform to the model out to a horizontal range R; at and beyond R, an increase or 
decrease in aerosol density exists. 

For all conditions the relationship between signal detectability, the altitude of the 
test vehicle, and the horizontal range R is given. 
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Several examples are provided to illustrate the use of the curves. Also presented 
are analyses of detection capabilities for the system used in verification of the model and 
a proposed improved system which is within the present state of the art. 


Derivation of Parametric Curves 

Condition I: total signa l det ectabil ity.- By using equation (3) the total number of 
photoelectrons produced at the photocathode of the receiver in counting time t c is 


N, 


e,t - 


^I ft (h)e- 2 W R 

4R 2 


( 12 ) 


where 


Q - 2A r t c yrj - 


A ryr ? 
AB 


The S/N value for total signal detection is 


S N e>t 

N |/N e>t + N e>D + N e>B 


(13) 


where N e g is the number of thermal electrons (dark current) and Ng^g is the num- 
ber of electrons due to sky brightness. It is assumed that » Ng^g + N e? g; this 

assumption is discussed in appendix F. Equation (13) then reduces to 



(14) 


Substitution of equation (12) into (14) yields 

s=f^ij e ft (h)R xio-i 

I R2 phj 


(15) 


Neglecting, for the moment, attenuation, both sides of equation (15) can be set equal to a 
parametric variable K: 





x 10" 14 


(16) 
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and 


K = 



(17) 


By using equation (16), K is plotted as a function of altitude for S/N = 1 and ^10 in 
figure 10. Figure 11 gives the relationship (eq. (17)) between QE-p and R for various 
constant values of K. By using the curves in figures 10 and 11, it is possible to 
(A) determine an existing system's detection capability, and (B) select system parameters 
to meet predetermined detection criteria. In determining a system's detection capabili- 
ties, a maximum range of the system is determined for a given QE T . This range R 
must be corrected for attenuation. Referring to equation (12), the relationship between 
predicted range R and range corrected for attenuation Rq is 

R = R c e^ h)R C 


This relationship is plotted for constant values of %(h) in figure 12. To determine R^ 
the value of /S^(h) from figure 13 corresponding to the value of h used in determining 
R is used. The intersection of this /3^(h) curve in figure 12 with the value of R then 
determines R^;. 

In selecting the system parameters to meet predetermined detection criteria, the 
quality factor QE-p of a laser system is determined for a given range R. This factor 
must also be corrected for attenuation effects. By referring to equation (12), the correc- 
tion to QE-p is 

(Q^t) c = QE T e 2(3t(h)R 

The correction factor e ^^t^ R j S plotted for constant values of 0t(h) as a function of 
range R in figure 14. /^(h) determined from figure 13 and the correction factor 
from figure 14 is applied to the predicted quality factor QE-p as above to obtain 

(qe t ) c . 

A procedure for making both of these determinations with examples is as follows: 

(A) Determination of an existing system's detection capability 

(1) Using known or measured parameters of the system, calculate QEt- 

(2) Select an altitude and minimum S/N value desired (S/N = 1 or \Jl0j, and 
refer to figure 10 to determine the corresponding value of K. 

(3) Refer to figure 11 and determine the value of R at which the constant K 
curve found in step (2) crosses the known value of QE-p. This value of R 
indicates the maximum detection range. 
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(4) Correct for attenuation as outlined above. 

Example 

(1) Given: QE-p = 6.61 x lO"-^ J-sec-m2 in the present system. 

(2) Select h = 0.31 km and S/N = |f!o. 

(3) By referring to figure 10, K = 1.6 x 10" 10 is obtained. By using this 
value in figure 11(c) (as illustrated), a maximum range of 5 x 10 2 m is 
obtained. 

(4) R c = 4.8 X 10 2 m. 

This example illustrates the low sensitivity of the present airborne laser system. 
The altitude in this example was chosen to match the altitude at which the sample oscillo- 
scope trace, figure 7, was taken. As can be seen in figure 7 at a range of 4.8 x 10 2 m, 
the S/N value is approximately ^10. 

(B) Selection of system parameters to meet predetermined detection criteria 

(1) Predetermined detection criteria specify altitude h, range R, and S/N 
value. 

(2) From figure 10, determine the value of K which corresponds to the pre- 
determined values of h and S/N. 

(3) From figure 11, determine the value of QE-p at which the line corre- 
sponding to value of K found in step (2) crosses the predetermined value 
of R. 

(4) Correct for attenuation as outlined previously. 

Example 

(1) Given: h = 10 km, R = 10^ m, S/N = /l0. 

(2) From figure 10, K = 4.8 x 10"10. 

(3) From figure 11, it is determined that the constant K = 4.8 x 10" 10 curve 
crosses R = 10^ m at a value of QE^ = 2.3 x 10"H J-sec-m2. 

(4) (QE t ) c = QE T (1.08) J-sec-m2 

= 2.5 x 10"H J-sec-m2 

Condition II: molecular signal detectability . - For this condition, it is assumed that 
there exists a clear atmosphere in front of the test vehicle to a range R, at which there 
is an instantaneous change in the molecular density with no change in the aerosol density. 
As in equation (12), the number of photoelectrons due solely to molecules can be written 
as 


£cQE T -2&(h)R 

N e ,R= ^ f R (h)e PtV 

4R 2 
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With a percent change in molecular density m, the signal-to-noise ratio is 


S m mN e,R x 10-2 

N ' fv 

This then gives 

— 6.17 

Q E T _ K = N » 12 

R 2 mf R( h ) 

neglecting attenuation. A family of curves for constant values relating m with altitude 
is shown in figure 15 for both S m /N = 1 and |flO. The procedure for utilizing these 
curves with an example is given below: 

(A) Determination of an existing system's detection capability 

The procedure is the same as that given for (A) in condition I, except for step (2). 
In step (2), select also a percent change in molecular density m and then refer 
to figure 15 for the corresponding value of K. 

(B) Selection of system parameters to meet predetermined detection criteria 

The procedure is the same as that given for (B) in condition I except in step (1) a 
percent change in molecular density m is also predetermined and in step (2) 
figure 15 is referred to for the corresponding value of K. 

Example 

(1) Given: h = 16 km, R = 10^ m, S m /N = \[l0, m = 1 percent. 

(2) From figure 15, K = 9 X 10“®. 

(3) From figure 11(a), it is determined that the constant K = 9 x 10"^ curve 
crosses R = 10^ m for a QE-p = 8.0 X 10“7 J-sec-m 2 . 

(4) With h = 10 km from figure 13 /3j-(h) = 2.0 X 10 - ® m _ l. By referring to 
figure 14, the correction factor for R = 10^ m is 1.08; therefore, 

(qe-p^ = 8.0 X 10-7(1.08) J-sec-m2 = 8.64 X 10“7 J-sec-m2 


Condition III: aerosol signal detectability .- This procedure is analogous to that of 
condition II. It is found that with attenuation neglected 



6.17^fc(hj 


x 10“ 12 
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A family of curves for constant values of K relating a with altitude is shown in fig- 
ure 16 for both S a /N = 1 and JlO. The procedures for utilizing these curves are given 
below: 

(A) Determination of an existing system’s detection capability 

The procedure is the same as given for (A) in condition I except a value for per- 
cent change in aerosol density a instead of m is selected and the value of K 
is obtained from figure 16. 

(B) Selection of system parameters to meet predetermined detection criteria 

The procedure is the same as given for (B) in condition II except a percent change 
in aerosol density a instead of m is predetermined and figure 16 is referred to 
for the corresponding value of K. 

It should be noted that although a distinction is made in conditions II and III between 
the molecular and aerosol density changes, the laser system measures only total scat- 
tering and cannot distinguish between the two. 

Use has been made of these curves to predict the maximum range Rq as a func- 
tion of altitude for two systems: the present system (QEqi = 6.61 x 10"^) and an 
improved system (QE-p = 3.75 x 10 -10 ]. Table III summarizes the characteristics of 
these two systems. Shown in figure 17 is the predicted maximum range Rq capability 
for signal detectability of a clear atmosphere for a S/N value of ]jld. In figure 18 the 
maximum range R^ as a function of altitude for the improved system is plotted for the 
detection of a 1-percent change in molecular density with a S m /N ratio value of 1; the 
present system could not make this detection for a range greater than 10^ meters at any 
altitude. 

Figure 19 gives the maximum range Rc for the improved system as a function of 
altitude for the detection of a 10-percent change in aerosol density with S a /N equal to 1. 
These design curves should be very valuable to investigators either presently involved in 
an airborne laser program or contemplating the initiation of such a program. They pro- 
vide a means of quickly determining either the capabilities of their present system or 
system parameters necessary to fulfill the objectives of their programs. 

CONCLUDING REMARKS 

A theoretical model of ruby-laser radiation backscatter with a wavelength of 
0.6943 pm in the lower atmosphere has been developed. This model includes not only 
molecular scattering based on Rayleigh theory but also atmospheric aerosol scattering 
based on Mie theory. 
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Based on this model and the measurements made, it appears that molecular scat- 
tering predominates in the altitude range of 4 to 14 km. 

This model has been verified up to an altitude of 11 km by using a laser system 
installed in a jet airplane. The system was found to be satisfactory for measuring the 
average scattering function and results were in good agreement with the model. 

The airborne system did not, however, have sufficient performance capability to 
probe the small horizontal scattering changes that might be associated with such phenom- 
ena as wind shear layers, temperature gradients, and possibly clear-air turbulence. A 
set of design parametric curves for an airborne laser probing system has been derived 
based on the theoretical model. These design curves should be very valuable to investi- 
gators either presently involved in an airborne laser program or contemplating the initia- 
tion of such a program. They provide a means of quickly determining either the capa- 
bilities of their present system or system parameters necessary to fulfill the objectives 
of their programs. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., September 12, 1969. 
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APPENDIX A 


DERIVATION OF GENERAL SCATTERING EQUATION 


The equation below gives a fundamental relationship between the power transmitted 
and the power received: 


p r _ CT R,( e > x ) N R( h ) + g M^> x )N M (h) ^ e -2/S t (h)R 
P T A t (R)R 2 


(Al) 


where R 2 » A r and Aq'(R). The total number of Rayleigh and Mie scattering centers 
can be related to Rayleigh and Mie number densities, respectively, by 


N R (h) = n R (h)v 

(A2) 

N M (h) = n M (h)v 

(A3) 


The scattering volume in terms of the laser pulse width is 


v = ctA-ji(R) 


(A4) 


The average power transmitted per pulse can be related to the total transmitted energy by 



(A5) 


Substitution of equations (A2), (A3), (A4), and (A5) into equation (Al) yields 


p r = 


cE-pAj 


2R^ 


°R( 0 >^)nR(h) + a M (0,X)n M (h) 


s -2 j 8 t (h)R 


where the factor of 2 in the denominator is due to the conversion of the time variable to 
oscilloscope time. 
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APPENDIX B 


DERIVATION OF RAYLEIGH SCATTERING FUNCTION 

Although the incident laser radiation is polarized, unpolarized scattering theory for 
8 = 0° or 180° predicts correct results; therefore, it will be assumed for purposes of 
calculation that the incident light is unpolarized. Thus, the cross section for scattering 
due to molecules in the notation of the present report can be shown in reference 13 to be 

a R (0, X)=A — — — x 10 24 m2sr-l (Bl) 

Ck 

The equation for the anisotropy factor from reference 13 is 

F _ 3(2 + D) 

6 - 7D 

where D has been found to be 0.03 (ref. 14); thus, F = 1.05. With 8 = 180°, 

A = 0.6943 pm, and a ~ 1.73 X 10 -30 n (ref. 15), 

ct r (180°, 0.6943 pm) = 2.10 X 10 -32 m2 S r~l (B2) 

Substituting equation (B2) into equation (2), the scattering function due to molecules 
is obtained as a function of molecular number density: 

f R (180°, 0.6943 pm, h) = 2.10 x IQ" 32 n R (h) sr -l m -l (B3) 
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APPENDIX C 


DERIVATION OF MIE SCATTERING FUNCTION 


Scattering due to aerosols is more difficult to analyze than that due to molecules. 
Aerosols consist of a variety of airborne particles and include water or fog droplets, 
dust, and combustion products. The difficulty arises because the size of these particles 
generally is on the order of the wavelength of the incident light, and resonances may 
occur. 

Again, by using the unpolarized incident-light theory, the cross section for scat- 
tering by aerosols is shown in reference 13 to be 


CT M( e > x ) 


i 1 (0,X,t,p) + i 2 (0,A,i,p) 
2k 2 


x 10“ 12 m^sr - ! 


(Cl) 


where ij and i 2 are the Mie scattering amplitudes and are a function of the angle 6 , 
the wavelength of the incident light A, the relative index of refraction of the particles t, 
and the shape of the aerosols p. The shape of the aerosols is assumed to be spherical. 
This assumption is reasonable since many of the particles are liquid, and these tend to 
form spheroids. It is also a necessary assumption since the scattering amplitudes have 
been calculated for spheres, but little work has been done for other shapes. Tabulation 
of Mie scattering amplitudes calculated before 1957 is given in reference 13. More 
recent tabulations for spherical particles are reported in reference 16. The average 
index of refraction will be assumed to be 1.50. (See, refs. 13 and 17.) This assumption 
is based on the discovery of dissolved salts in the liquid aerosols. (See ref. 18.) Based 
on these assumptions and with 6 = 180° and A = 0.6943 /um the cross section for scat- 
tering due to aerosols is 


ct m ( 180°, 0.6943 (im) = 


4(180°, 0.6943 (im, 1.50, r) + i 2 (180°, 0.6943 fxm, 1.50, r) 

2k 2 


X 10"12 m2sr~l 


(c 2) 


In order to use equation (C2) in equation (2) of the text, it is necessary to know the 
size distributions of the atmospheric aerosols. Investigations of the size distribution of 
atmospheric aerosols have been performed by Junge (ref. 19) and others (see, for exam- 
ple, refs. 20 and 21). Junge's measurements were made over the European Continent but 
agree well with data taken near Baltimore, Maryland (ref. 20). The distribution from 
r ~ 20 pm to r ~ 0.1 jam is given as 


dn M (r,h) = C(h)r -3 d log r 


(C3a) 
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APPENDIX C 


then 


d log r = 0.434r~l dr 

and 

dn M (r,h) = C'(h)r-4 dr (C3b) 

This distribution has been reported by Junge (ref. 22) to be invariant with altitude in the 
lower atmosphere. 

Substituting equations (C2) and (C3b) into equation (2) and integrating yields 


f M (180°, 0.6943 pm, h) = 


J r-, 


21r 


(C4) 


By changing the variable from the radius r to a nondimensional size parameter g, 
where g = 27 jt/X and substituting k = 2 tt/A, it can be shown from equation (C4) that 


f M (180°, 0.6943 pm, h) = C’(h)| f 

A %) c 


g 2 ijtg, 1.50) + i 2 (g, 1-50) 


dg x 10“ 6 sr-lm~l 


Si 


= C'(h)^I x 10" 6 
A 


(C 5) 


where 


I = 


§2 ii(g» L5 °) + 1-50) 

§1 g 4 


J. 


dg 


Integrals similar to that in equation (C5) have been evaluated numerically in ref- 
erence 17 for a number of wavelengths. These results over a size range of r j = 0.04 pm 
to r 2 = 10 pm are given to the third significant place as 


I = 0.442 (C6) 

The constant C T in equation (C5) can be evaluated in terms of the total aerosol 
number density by integrating equation (C3b) over the size limits consistent with data 
obtained near Chesapeake Bay (ref. 23): 
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or 


n20 p20 

n M^) = \ dn M (r,h) = C’(h)\ r-4 dr (C7a) 

J 0.1 J 0.1 


C'(h) = 3n M (h) x IQ' 3 


(C7b) 


By substituting this value for C’(h) and the value of the integral from equation (C6) into 
equation (C5), the scattering function due to aerosols becomes 


f|^(180 o , 0.6943 jum, h) = 6.0 x 10~^ n^(h) sr'^m"! 


(C8) 


22 



APPENDIX D 


AEROSOL NUMBER DENSITY AS A FUNCTION OF ALTITUDE 

In 1961, a sea-level measurement near Chesapeake Bay of the total aerosol number 
density was made. (See ref. 23.) This measurement was performed by fitting light- 
attenuation data to theoretical predictions. The theory was based on the same assump- 
tions as presented in appendix C, that is, spherical aerosols with index of refraction of 
1.50 and Junge's size distribution from rj = 0.1 pm to ^ = 20 pm. The results taken 
on a clear day with a visibility of 25 km yielded a total aerosol number density of approx- 
imately 215 particles/cm3. 

Total number density as a function of altitude has been measured. All measure- 
ments made before 1954 are summarized in reference 24. It was found that the density 
decreased exponentially with altitude up to 5 km. The scale height of the decrease in 
aerosols was determined to be 1.2 km. Aerosol densities up to an altitude of 26 km were 
measured in 1965. (See ref. 25.) The aerosol number density at kilometer intervals is 
shown in table I. The aerosol number density from reference 25 was corrected since 
aerosols with radii greater than 0.125 pm were measured; and to be consistent with the 
densities from reference 24, it was necessary to correct these values to a lower limit of 
0.10 pm by assuming the Junge distribution given in equations (C3a) and (C3b). Because 
the value for the number density at 5 km differed in the two investigations (refs. 24 
and 25), the average value at this altitude was used. 

Another tabulation of aerosol number density as a function of altitude is presented 
in reference 26. The aerosol data which were used herein were obtained from refer- 
ence 27 for an altitude range of 10 to 25 km. These data are an order of magnitude below 
those used in the present report from reference 25. The justification for using the data 
of reference 25 is that it predicts an increase in scattering greater than Rayleigh predic- 
tion above an altitude of 16 km; this result has been observed by a number of investiga- 
tors. (See refs. 6, 28, and 29.) Data plotted in figure 20 from reference 29 is typical of 
the increased scattering in the 16- to 26-km altitude region reported by the other two 
investigators (refs. 6 and 28). 
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APPENDIX E 


DETERMINATION OF ATMOSPHERIC EXTINCTION COEFFICIENT 


The total atmospheric extinction coefficient can be written as 

0t(h) = %(h) + jS M (h) + /3 0 (h) (El) 


j8j^(h) can be calculated by integrating the Rayleigh scattering function over all space as 



decree, A) sin 0n R (h) 


(E2) 


or by combining equations (Bl) and (E2) 

fe(M - I 2 ’ ™ C ° S ^ )k4 ' 5 | 2 F„ R (h) X 1024 

which, for A = 0.6943 / am, gives 


(E3) 


/3 R (h) = 17.6n R (h) x 10“ 32 m _1 (E4) 

Now, at sea level /3q(0) ~0 and /3(-(0) -0.135 x 10“3 m - l (from ref. 23 for a visibility 
of 25 km); therefore, from equation (El) 

/3 m (0) = 0.131 X 10" 3 m-1 (E5) 


since 


/^(O) = 4.24 X 10“ 6 m" 1 

An effective extinction cross section for aerosols can be calculated from 

/3m(°) = n M(°) CT M 


(E6) 


(E7) 


where cr^ is the Mie extinction cross section. Based on data from table I, where 
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°M = 6.09 x 10" 7 cm 3 m - l 

The extinction cross section as a function of altitude then is 

/3 M (h) = ,6.09n M (h) x 10“ 7 m'l (E8) 

and the total extinction coefficient obtained by combining equations (El), (E4), and (E8) is 

/3 t (h) = 17.6n R (h) x 10 -32 + 6.09n M (h) x 10' 7 + j3 0 (h) m~l 

where n R (h) is obtained from reference 12, n M (h) is obtained from table I, and /3o(h) 
is obtained from reference 25. 


25 


APPENDIX F 


ON THE ASSUMPTION N 6j t » N e> D + N e> B 


The number of electrons at the photocathode due to dark current in a counting time 


is 


N* 


«J"a D 

q = — L~ t c (6.25 x 1018) electrons 
5 G 


or 


N P d = —^^(6.25 X 10 l 8 ) electrons 
e ’ u 2BG 


The number of photoelectrons produced at the photocathode of the receiver by the daylight 
sky brightness is 


N e , B = bfi r Ai^A r t c ? 7 yC X 10.76 photoelectrons 


or 


N e b = bOpA^Qj^ x 5.38 photoelectrons 

where b is the sky brightness in W-ft" 2 _sr _ l-A“^ and is equal to 10 "^ in the 6000- to 
7000- A region (ref. 30). The additional parameters for the present and improved systems 
are shown in table III. 

For the present system, Ng^g + N 6; b = 4.95 x 10"^ electrons and for the improved 
system, N e g + N 6j b = 8.88 x 10“2 electrons. For a signal-to-noise ratio of 1 or \/l0, 
N 6j t = 1 or 3.16 electrons, respectively. Therefore, for both systems and both signal-to- 
noise ratios the inequality N e> t » N e> g + Ng^g is a valid assumption. 
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TABLE H.- CHARACTERISTICS OF THE TRANSMITTER AND RECEIVER SUBSYSTEMS 


Transmitter subsystem: 

Laser transmitter Ruby 

Wavelength, jiim 0.6943 

Peak power output, MW 10 

Total energy output, J 0.1 

Q-switching technique Saturable filter (cryptocyanine in nitrobenzene) 

Output pulse duration, nsec 20 

Beam divergence after collimating optics (total), mrad 1 

Pulse repetition rate, pulse/sec 1 

Output monitor Photodiode 


Receiver subsystem: 

Receiving optics: 

Effective aperture, mm 110 

Acceptance half-angle, mrad 3 

Interference filter: 

Center wavelength, fim 0.6944375 

Optical bandwidth, pm 0.001175 

Transmission, percent 15 

Photomultiplier : 

Overall voltage, volts 1800 

Stages 16 

Quantum efficiency ^0.003 

Photocathode response S-20 

Photocathode diameter, mm 4.83 

Additional amplifier gain 10 
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TABLE HI.- PARAMETERS OF PRESENT AND IMPROVED SYSTEMS 


Parameter 

Present system 

Improved system 

y 

0.15 

0.60 

V 

2.8 x 10“ 3 

2.5 x 10-2 

A r , m2 

9.48 x 10-3 

0.125 

B, sec - l 

50 x 10 6 

5 x 10 6 

Q, sec-m2 

7.97 x 10~ 14 

3.75 x 10-10 

E-p, joules 

0.083 

1 

QE-p, J-sec-m2 

6.61 X 10-!5 

3.75 x 10-10 

JA Dj amperes 

6 x 10“ H 

6 x 10-H 

G 

8 x 10 5 

8 x 10 5 

Total field of view, mrad 

6 

0.2 

fi r , sr 

28 x 10" 6 

6.27 x 10" 8 

Av, A 

11.75 

2.0 
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Figure 3.- Block diagram of the airborne laser system used in this investigation. 
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Figure 4.- The laser system installed in 
(T) transmitter and receiver subsystei 





(T) Laser transmitter (4) Output monitor 

( 2 ) Storage capacitors (5) Transmitting optics 

(T) Laser power supply ( 6 ) Receiver 

( 7 ) Receiver power supply 

Figure 5- Transmitter and receiver subsystems. L-67-6365.1 
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Range, R, m 


® 


Peak of the integrated output of the photodiode, which indicates total energy 
output of the laser 

Crossover of the transmitted beam and the field of view of the receiver 
_2 

The normal R falloff of the atmospheric return. 


Figure 7.- Oscilloscope trace taken with airplane flying in level flight at an altitude of 0.31 kilometer over 
Williamsburg, Virginia (lat. 37°17' N, long. 76°45' W). 
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Total scattering function, f f (h), sr' 



Figure 8.- Comparison of experimental results obtained on the night of March 9, 1967, over Williamsburg, Virginia 
flat, 37°17' N, long. 76°45' W), with the theoretical model. 
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Correction to QE j, (QEj)^, , e^Ptt^JR 
Figure 14.- Attenuation correction to QEj for a given R and ftth). 
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Maximum corrected range, R^, m 


Figure 17.- A comparison of the present and improved system's limit of signal detectability with S/N - tjlo. 
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